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PREFACE 

The David W. Taylor Lectures were conceived to honor our founder in 
recognition of his many contributions to naval architecture and naval hydro- 
dynamics. Admiral Taylor was a pioneer in the use of hydrodynamic theory 
and mathematics for the solution of naval problems. He established a tra- 
dition of applied scientific research at the "Model Basin" which has been 
carefully nurtured through the decades and which we treasure and maintain 
today. It is in this spirit that we have invited Prof. Louis Landweber to 
be a David W. Taylor Lecturer. 

Prof. Landweber was born in New York in 1912. He received his Bache- 
lor of Science degree in Physics from the City College of New York in 1932. 
That year he came to the U.S. Experimental Model Basin at the Washington 
Navy Yard as a junior physicist. When the David Taylor Model Basin opened 
at Carderock in 1940, Dr. Landweber headed a small research group which 
shortly afterwards was expanded into the Hydrodynamics, Division. During 
his years at Carderock, he continued to head this division, and he also 
completed work on his Ph.D. at the University of Maryland. He was an advo- 
cate of advanced training for the staff and taught courses for the 
University at the Center. 

After twenty-two years of distinguished research at the "Model Basin," 
Dr. Landweber left to become a professor at the University of Iowa. He 
continued his research in ship hydrodynamics at the school's Institute for 
Hydraulic Research. During the intervening years since leaving the Center 
in 1954, he has maintained close ties with his colleagues here at the 


Center and has returned frequently for meetings and panel sessions. 
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FOREWORD 

After a lapse of twenty-four years, it was surprisingly easy and 
pleasant to readjust for a month to a daily schedule of work in a mezzanine 
office at the David W. Taylor Naval Ship Research and Development Center. 
I found Ship Hydrodynamics to be alive and well, its problems being vigor- 
ously attacked by a dedicated and talented staff, as in the "old days." 

The adjustment was probably equally easy and pleasant for another 
"alumnus" and David Taylor Lecturer, John Wehausen. It is remarkable that 
so many of the alumni have stayed in the field of Ship Hydrodynamics, an 
indication of the challenging nature and attractiveness of the subject. 

To all who visited me, to consult, gripe, reminisce or simply to 
educate me, who extended generous hospitality in accordance with the 
fortune-cookie admonition, "Take time to play in order to have a long 
life," who invited me and arranged that day-to-day living would be con- 
venient and comfortable for my wife and myself, in other words, to all 
who helped disprove the old adage that "one cannot go home," our deepest 


appreciation. 
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ON IRROTATIONAL FLOWS EQUIVALENT TO THE 
BOUNDARY LAYER AND WAKE 


INTRODUCTION 

One of the simplest and most useful results of boundary-layer theory 
is that the flow exterior to the boundary layer, which will be assumed to 
be irrotational, is pushed outwards by an amount called the displacement 
thickness. This suggested, early in the development of boundary-layer 
theory, that the accuracy of a boundary-layer calculation for a body could 
be improved by adding the boundary-layer thickness to the body dimensions 
and using the predicted pressure distribution on the so-thickened body in 
a recalculation of the boundary layer. The irrotational field about the 
thickened body, including the displacement thickness of the wake, is itself 
of great interest, and numerous attempts have been made in this manner to 
calculate the effect of the boundary layer and wake on the outer irrota- 
tional flow. 

The concept of the thickened body gives an approximate model which is 
usually justified by its consistency with the approximations of thin 
boundary-layer theory. We shall review the basis of this model and suggest 
ways of refining it. Such a development would be useful for several current 
problems of ship hydrodynamics, among them the determination of Betz 
sources in a method of calculating viscous drag by an analysis of wake 
survey data, and the investigation of the effect of the boundary layer and 
wake on wavemaking resistance. 

We shall assume that the boundary layer and wake (BLW) are known, 
and seek an irrotational model which yields the actual outer flow. An 
approximate solution of this problem has been given by Preston [1] and 


Lighthill [2]. 


Solution of Preston and Lighthill - Two-Dimensional Case 


Let us consider a two-dimensional flow about a body of small curvature 
in a uniform stream of velocity Us» on which a thin boundary layer is 
present. Coordinates parallel and normal to the surface will be denoted 


by (s,n), and the corresponding velocity components by (u,v). We shall 


suppose that the radii of curvature are so large compared with the boundary- 


layer thickness 6, that the curvilinearity of the coordinate system may be 


ignored. The normal velocity component at n = 6 is then given by 
6 6 6 
ve) = | Ban- [ Ba- 6 ‘| (U-u) dn -6 
0 0 
6 
yah rst ow (-¥)a (1) 
as See a uy) oe 
0 


where U = u (s,6) and 7 is the displacement thickness. In the irrotational 
flow about the body, with velocity components U,(s.n), V,(s.n), we have 
(0U,/dn) 4 = 0 where the subscript 0 denotes n = 0. Hence we may put 


U,(s,n) = U,(s), since 


dU 
o pee leas - 
U,(s,n) = U, (s,0) +n (2 I, 5 U,(s,0) = U,(s), O< a <6 


The normal velocity in this irrotational flow, at n = 6, would then be 


given by V,(s,9) = V,(s,0) + 6(8V,/0n) 9, and hence, since V,(s,0) = 0 and 
dV, /on == dU, /9s, we obtain 


dU, 
V,(s56) = - 6 Ace (2) 


We also assume that U = u(s,4) = U,(s,9) for a thin boundary layer. Thus 
the first term of the right member of (1) is attributable to this irrota- 


tional flow, and the second term represents an additional outward flow due 


to the boundary layer. 


Rute WY S <- (US) and consider that the "known" values of V"' on the 


contour n = 6, the edge of BLW, pose an exterior Neumann boundary-value 
problem for determining a source distribution m(s) on the surface of the 


body and (if the body is symmetric) along the centerline of the wake. 


Preston and Lighthill give the solution without hesitation. Their result is 


plese _i a 
eS oa (s,6) = on de (U5) (3) 


for which the argument given by Lighthill is: "This additional outflow is 
exactly 'as if' the irrotational flow around the body were supplemented by 
the effect of a surface distribution of sources, whose strength ... is ..." 
Equation (3) does not seem as obvious to the present writer. In order 
to derive it, one must assume that the velocity U(s) is undisturbed by the 
distribution m(s), an assumption which may not be consistent with the 
uniquely defined solution of the Neumann problem. According to this assump-— 


tion, the body surface is an equipotential for the distribution m(s) (be- 


cause the tangential velocity due to m(s) is zero). This immediately yields 


+ 
m(s) = + v'(s,0,) PE (4) 
z 277 Sober 0 
v' = 0 
since V'(s,0_) = O on the interior side of the equipotential surface and 


v'(s,0°) - V'(s,0 ) = 2mm(s). Then, from the Taylor expression V'(s,6) = 
v' (s,0,) + 6 (=~) 0 +... = V'(s,0,), we obtain the approximate solu- 
fenton ((3})) 4 
The special case of a flat plate, (even in a nonzero pressure 
gradient), however, confirms the approximate solution (3) without requiring 
the equipotential assumption. Since the plate is the limiting case of a 
y! 
m(s) 
m(s) 


Figure 1 


very thin body, the distributions m(s) are present on both sides. Applica- 
tion of the Gauss flux theorem (the usual proof) immediately gives Equa-— 


ieaoyn) (8) 4 


Another case that yields (3) is that of a constant source distribution 
on a circle. This gives a velocity potential due to a source of strength 


21™ma at the center of the circle of radius a, and hence V'(s,n) = V'(s,0) = 

(21ma) 
a 

to V'. Hence, if we consider a nonuniform distribution which is large in 


= 21™. But the local element about the point s contributes Tm 


the neighborhood of s, and negligible elsewhere, we would obtain V'(s,n) 5 
TM, half of the Preston-Lighthill formula. This suggests that, without the 
equipotential assumption, Equation (3) would give a good approximation only 
for very thin two-dimensional forms. 

The distribution m(x) on the surface of the body y = f(x), which 
satisfies the boundary condition do/ dn. = v'(s,0) can be determined as the 
solution of the Fredholm integral equation of the second kind, (see 


Figure 2) 


Tm(s) + § m(t) = in 18 dt = V'(s) (5) 
s 


in which t also denotes arc length along the body contour. The indetermi- 


nacy of the kernel 2. Qn Te at s = t can be removed by applying, the 


identity (at a smooth point of the contour) 


p) _  @ ri 
f sine, at = fF o,, ae =o (6) 


where the complex vector 2g e » extending from s to t, has given the 


pair of conjugate functions 


to which the Cauchy-Riemann Equations have been applied. We then obtain 


from (5) and (6) 


21™m(s) + ' [aco 52 an Top 7m(s) a Ln se | dt = V'(s) (7) 
s t 


Since the new integrand vanishes at the point s = t, where the largest con- 
tribution to the integral in (5) occurs, the form in (7) suggests that the 
Preston-Lighthill formula should give a good first approximation. 

From formula (3) for the source distribution, Lighthill immediately 
concludes that the displaced surface, obtained by the addition of the dis- 
placement thickness, is a stream surface. This again requires the assump- 
tion that all the flux from the source distribution m(s) be outward, i.e., 
that the given contour be an equipotential for the distribution. For then, 


the Gauss flux theorem gives the relation 


s 
an | m ds = U n(s) 
0 
where n(s) denotes the stream surface generated by m(s). Comparison with 


(3) then shows that n(s) = 61: 


A Second Approximation for the Source Distribution —- 
Two-Dimensional Case 

It appears to be easier to derive a formula for a source distribution 
M(x), equivalent to the BLW, for a thin, symmetric body, when M(x) is dis- 
tributed on the axis of symmetry, the x-axis. The given profile will be 
denoted by y = f(x), the edge of the boundary layer, EBLW, by y = g(x), and 
the profile thickened by the displacement thickness by y = h(x). Put 


Figure 3 


df/dx = tan y. For the irrotational flow, the velocity components in the 
x- and y-directions will be denoted by u(x,y) and v(x,y). 
We wish to determine the source distribution M(x) such that, on 


y = g(x), the stream function (x,y) assumes the values 


p = U(6-d,) (8) 
We have 
But 
WAGES) = WIG%5O) Pe ae) = T(x) - g u Gx, 0) (10) 
BUG) Se wC<5,O)) ar ws Q) = u(x, 0) (11) 


since a) = 0 by symmetry. Also, by (8), we have 


d dx 
= WG-8 ae Ne (12) 


Ss 


a 

Q)}Q 

n 

Se 
i] 


and hence, substituting (10), (11), and (12) into (9) and solving for M(x), 


we obtain 
m™(x) = $= [g(x) u(x, 0)-U(5-8,)] (13) 


Since, as is seen from Figure 3, 
g(x) = f(x) + 6 sec y (14) 
and 
U = u(x,g) sec y = u(x,0) sec y (15) 


we can write (13) in the form 
T(x) = = {u(x, 0) [£(x)+6, sec y]} (16) 


Here we have assumed that f(x), 6, 6 (u-U5); v and y are small quantities 


iL? 
of the first order and terms of third and higher order have been neglected. 

For the irrotational flow about the body without a boundary layer, (16) 
yields the well-known, second-order approximation for a centerplane 


distribution, 
m™ (x) = S [u(x,0) £(x)] (17) 
0 ~ Gb 2: 
Hence, the additional source strength due to the boundary layer is given by 
MG MG) = = = suGs0)) 6H see yl (18) 
0 “17 dx % iL 


The .result in (18) resembles the Preston-Lighthill formula (3), especially 


when the body slope is small. The distribution strength is approximately 


doubled because a pair of distributions on the body surface has coalesced 
on the centerplane, as was done in the example of the boundary layer on a 
flat plate. 

Similarly to (14), the displacement-thickness profile, y = h(x), may 


be expressed by 


h(x) = £(«) + Oy sec Y (19) 
and (17) becomes 
MG) 5 = & [u(x,0) h(x] (20) 


Comparison with the form of (16) now shows that, to the second order of 
accuracy, the irrotational flow about the displaced surface satisfies the 


specified boundary condition (8) at the edge of the boundary layer and wake. 


Solutions for Axisymmetric Flow - First Approximation 


Let TQ (s) denote the radius of a body of revolution, where s is arc 


length along a meridian section of the body, and r the radial distance to an 
arbitrary point. Let n denote distance from the body along the outward 
normal to its surface. We shall also employ cylindrical coordinates, (x,r) 
where the axis of symmetry is taken as the x-axis, and put TQ (s) = £(x) 


and tan y = df/dx. The edge of the boundary layer is defined by the surface 


& (x) 
n 
s £(x) 
6' 
T 
x = Se SS OY Sim 47% 
=O PO 
ore 
Figure 4 


ro(s) = ry (s) + § COS 7 Or BG) = FG!) + Oo GY) cos vy G)) (21) 


By a Taylor expansion about x, with dy/ds = - k, where k is the curvature 


of the meridian profile, and with x' - x = 6(x%') cos y(x'), g(x) becomes 
f ees 
g(x) = £(x) + 6 sec y + 66' tan y - x ke (eaioy 7 SEC W4P coo (22) 


The element of arc along a curve n = constant is hy ds = (1+kn) ds. The 


equation of continuity is then 


0 (ru) 


fe) 
Se ——— 45 = 
a + a [ (1+kn) rv] 0 (23) 
where u(s,n) and v(s,n) are the actual velocity components in the directions 
of increasing values of s and n. We shall also require the velocity com— 
ponents U(s,n) and V(s,n) of the "equivalent" irrotational flow. 
The boundary-layer thickness of an axisymmetric boundary layer is 


usually defined as 
2 eo po B@om) 
6, [2 E Te |e (24) 
0) 
An alternative definition, in which a higher-order term is retained, is 


given by 


6x 6 
r U(s,6) dn = J r[U(s,6)-u(s,n)] dn = ro U Oy 
0) 


Since r=rfr 


0 +n cos y, this yields the quadratic equation 


ry 6* + 6x? cos Y = rp 6, 3 (25) 


and the relation 


§* = ——____________. (26) 


The equivalent irrotational flow is determined by the values of v(s,é). 


By (23), we have 


6 
(14+k6) re v(s,6) = - if = [xr (U-u)-rU] dn 
0 


| 
S——),. 
QQ 

| 
= 

i] 

[<j 
Ne 
Qu 

t=} 

I 


if 2 
dean gece d Omed 
Sas Sole 8 ae SY) = arg eS 


where U = U(s,6). In the irrotational flow about the body without a bound- 


ary layer, with velocity components (U,5V})> (27) becomes 


6 


i $— (ru, (s.n)] dn 
0 


(14+k6) rs V,(s,9) = 


2 
ss d oO. al 
=- 6 AE (r Uy) Ftoetas (U, cos y) 


: ze oF Dixtarah 
since (0U,/dn) 9 = 0, and U,(s,n) = U, (s,0) +0(6°) = U,(s). If the 
boundary layer is thin, then U,(s) = U and the first term of the right 
member of Equation (27) can be interpreted as that contributing an addi- 
tional flux due to the boundary layer. We define the additional normal 


velocity V' due to the boundary layer by 


(14k6) rv" = (x )u6 (28) 


D 


10 


If the body were a circular cylinder of constant radius o> we could 
immediately determine the source strength m(s) equivalent to the effect of 


the boundary layer. Application of the Gauss flux theorem gives 


! => 
271 .V 4m (21x 9m) 


and hence by (28), 


eal el 
mS Te aig Sad (29) 


A generalization of (29) can be derived from the integral-equation 


formulation of the Neumann problem for the prescribed boundary-layer flux. 


If P and Q are points on the body surface, we have 


2mm(P) = il m(Q) s-s+ dSy + V'(P) (30) 
5 Q “PQ 


The singularity of the kernel when Q coincides with P can be removed by 


means of the relation, valid at smooth points of the surface, 


|e == eg (31) 


Thus (30) may be written in the form 


pil 


4rm(P) = i [aco se - me) 5 


ae |e + V'(s,0) (32) 
S Ee “Q * 


which suggests the first approximation 


iLL 


1 


m(P) = = WSs Ores 


(x 9U5,) (33) 


a generalization of (29). 

In order to improve upon the approximate source distribution (33), we 
shall now consider an axial source distribution M(x). For irrotational 
flow about the given body, r = f(x), the modified Munk formula of Landweber 


[3] gives the approximate solution, in terms of the free-stream velocity U 


0’? 
My) 5% Gtk.) Uy & (34) 
1g 0 dx 
where ky is the longitudinal added-mass coefficient. For the thickened 
body, with r = r, (x) = £(x) + O5 sec y, formula (34) gives 
1 d 2 
M, (x) BiG (1+k,) U5 dx ‘1 (x) 
=M +5 (+k, ) U (46 sec Y a5 6 2 Bee” Y) (35) 
~ 0) = 1 2a: 


If the alternative displacement thickness 6* is used, the same expression 


is applicable, with 65 replaced by 6%. 


Axisymmetric Flow - Second Approximation 


We shall now derive a second-order approximation for an axial source 
distribution for slender bodies. The given profile is r = TQ (x) = 35 (6%), 
the edge of the boundary layer is r = g(x), and the displacement-thickness 
profile is r = h(x). The velocity components in the x- and r-directions 


for the equivalent irrotational model are u(x,r), v(x,r). 


12 


Figure 5 


Let us apply the Gauss flux theorem to the volume of the body of 
revolution generated by the profile r = £(x), extending from the nose to a 
transverse section at x, assuming that the flow about the body is generated 
by an axial source distribution My (x) - Since r = f(x) is a stream surface, 


the flux theorem gives 


2mru, (x, 4) dr (36) 


a 

=| 
co 

Ss 

i=) 

a. 

tes 

Il 
oo 


Here [u, (xr), v,(s4)] denote the velocity components of the irrotational 
flow when the given profile is a stream surface, as distinct from the 
components (u,v) of the equivalent irrotational flow. Differentiating (36) 


and integrating by parts, we then obtain 
if 


IL al 2 iL 
My Cs) = Te ax | fs u, Gx, £)- il ic jam ar | Gi) 
0 
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We shall assume that f, 6 and y are small of first order, and neglect 


du, d(rv,) 
terms of fourth and higher order in (37). Then, since r Been ae 
dM 
= 2r = is small of order a (37) yields 
dik dade 
My (x) are he Giz Lie u, Gx, £) ] 
or, since u, Gs, £) = U, (s,0) cos Y = U,(s,5) cos Y 
1d 2 
My (x) ite [£ U,(s,9) cos y] (38) 


Similarly, the source distribution M, (x) for the irrotational flow about 


the surface displaced by the displacement thickness 6* is given by 
M, (x) = & [ (£+6* sec y)- U_*(s,6) cos y] 
1 4 dx i ? 


where U,* denotes the s-component of the velocity. Then, assuming that 


U,(s,9) = U,*(s,6), we obtain 


M, (x) - My Gd) = + 2 [(2£6*+6%" sec y) U,(s,6)] (39) 


Similarly, when the surface of revolution is generated by r = g(x) and 


the axial distribution is M(x), the Gauss flux theorem gives 
x s g 
AT [ M(x) = i 27g(x) V(s,6) (1+kd) dS + 27ru(x,r) dr 


0 0 0 


or, by (27), and integration-by-parts of the last integral, 


14 


2 
© a d 2 
Mic) = a (UE5,) - 6 a (uf) - $- & W cos y) +$% ioe 1 


The same analysis, applied to the source distribution My x) » yields 


2 
Py d Om ids 
2M (x) Oe ae 6 dx (U,f) ae D dx (U cos ) wD) dx lg AE (x, g)] 


IL 
and hence, assuming U(s,6) = U,(s,6) and u(x,g) = u, (x, 8), we obtain 


M(x) - My (x) = $a (UES, ) (40) 


Comparison of (39) and (40) now shows that M, (x) = M(x) if 


£6% + Sox sec y = £6 (41) 


iL 


This agrees with the definition of 6* in (25) when y is small. Otherwise, 
the displaced surface should be taken in accordance with (41) rather than 
(25). 


A Vorticity Theorem 

Consider a steady mean flow of an incompressible fluid about a body at 
rest in an otherwise unbounded fluid, with mean vorticity present in the 
boundary layer and wake, BLW, bounded internally by the body surface S and 


externally by the surface T; see Figure 6. The boundary conditions to be 


T 
U : 
————_ > 
+ BLW 


Figure 6 
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satisfied are: 

1. The velocity vector von S45 the exterior side of S, satisfies the 
nonslip condition, v = 0. 

2. The flow exterior to T is irrotational. 

Then we have the following theorem: 

The disturbance flow field exterior to S can be generated by the vorticity 
in BLW alone. 

Proof: 

First let us suppose that distributions of vorticity and irrotational 
singularities are present in B, the interior of the region bounded by S. 
The velocity components normal to S induced by these distributions define 
a Neumann problem for a source distribution M on S, equivalent to the in-— 
terior distributions. 

We now have a source distribution M on S and vorticity W in BLW. The 
tangential velocity Wes is continuous across S, and hence, by the nonslip 
condition, Vis 0 on the interior side of S. Also the flow induced by the 
distributions M and W is irrotational within B. Hence v = 0 within B, and, 
consequently, the normal component of the velocity on the interior side of 
S is zero. But, because of the nonslip condition, va is zero on the ex- 
terior side of S also, hence the strength of the source distribution must 
be zero. This leaves only the vorticity distribution W, as we wished to 


show. 


Verification with Stokes Solution for a Sphere 
Stokes solution for a sphere of radius a in a uniform stream in the 


x-direction, expressed in spherical coordinates (R,9,>), is given by 


Figure 7 
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PO |b 


v, =U cos 8 + 2 (er 
Ro 


) cos 6 (44) 


v, = - U sin 0 + ae sin 0 (45) 
8 R? R 


Here the first and second terms represent components due to the uniform 
stream and a doublet at the center of the sphere, and the last terms are 
due to vorticity. We shall now show that the velocity at points of the 
x-axis, (8 = 0, R > a), can be obtained from the vorticity outside the 
sphere, together with the uniform stream. 

The vorticity, determined from the last terms of (44) and (45), is 


given by 


Wp = Wy = 0, Ww, = 28 shim IS) (46) 


At constant R and 8, a vortex ring of unit strength induces a velocity (by 


the Biot-Savart Law), at a point x of the x-axis, given by 


R- Sac ) 
Neem Gee Gr 


The velocity due to the vorticity in the space exterior to the sphere is 


then given by 


eee J i oa sin 9 dé aye (47) 
4% (x -2xR cos 6+R ) 
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With the substitution U = cos 9, integration with respect to wu yields 


u=-—4, i 4 joea’? aN z ) s2can?) [n-x| tox) 
an |R-x| | R+x| 


i 
4 ([Rox|) "= (Re) * | dR 


This is readily evaluated by considering separately the integrations from 


a to x, and from x to ©. The former yields 


x 
ANB _ _ Wa Dad 
ie een | RdR = mares (x -a’) 
x 2x 
a 
and the latter gives 
u-fa f Bo 
R 
x 
Hence 
= _ Ua 3Ua 
u = uy + uy Ti aaghs oS (48) 
2x 


in agreement with the distrubance velocity in (1) when 9 = 0, given in (44). 
This result indicates that the doublet at the center of the sphere is 
an irrotational equivalent of the negative of the internal vorticity, since 


the sum of their fields is zero exterior to the sphere. 
Some Applications of the Vorticity Theorem 


We have previously investigated irrotational equivalents of the bound—- 


ary layer and wake such that the outer irrotational flow was preserved. It 
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was seen that this could be accomplished in various ways, such as a source 
distribution on the given surface and along the wake, a source distribution 
on the edge of the boundary layer and wake, or by the irrotational flow 
about the displacement-thickness surface. In none of these irrotational 
models was it possible to retain the given body as a stream surface. 

Intuitively, it appears desirable to match the boundary conditions on 
both the body and the edge of the boundary layer and wake, in order to 
obtain a more realistic irrotational model. We have seen that the vorticity 
distribution alone yields both the nonslip condition on the body and the 
boundary condition on EBLW. In the previous models only the boundary condi- 
tion on EBLW was employed. An additional source distribution, on the body 
or in its interior is required in order to satisfy the condition that the 
body remain a stream surface. The condition of zero tangential velocity 
would not be satisfied, but this seems to be physically less important in 
an irrotational model. 

The boundary conditions on the body surface S and on EBLW define a 
Neumann problem which can be readily formulated as a pair or integral equa- 
tions to be solved simultaneously for a pair of source distributions. The 
locations of the source distributions may vary, even for a given body, as 
has already been illustrated. If these are taken to be distributions m(P) 


on S and U(Q) on T; one obtains the integral equations 


0 ill ) il 
27™m(P) — [ac oe dG = fra Gor dS,, 
on, Tppr P on, TpQ! Q 
S T 
es Ox 
ms Uo dn, 9) 
e) 1 ) 1 
27U(Q) - [va ea éSag = [ sen aa = By S WQ) Go) 
} ong ToQ! Q ; ony T p19 P 


where P and P' denote points on S, and Q and Q' on EBLW, and v(Q) denotes 


the normal component of the velocity at T, which is assumed to be known. 
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The Poincaré Transformation [4] 

It is well known that the velocity field induced by a vortex ring is 
identical to that of a doublet sheet on a surface capping the ring. The 
Poincaré transformation, derived below, enables more general relationships 
between fields of vorticity and irrotational distributions to be obtained. 
Since the entire disturbance flow field is induced by the vorticity alone, 
the velocity field of the equivalent irrotational distributions would be 
identical to that induced by the vorticity in the regions exterior to the 
vorticity domain. For the case where the vorticity lies in the boundary 
layer and wake of a body, the vorticity induces not only an irrotational 
field exterior to its outer boundary T, but also an irrotational field 
within the body. Because of the nonslip condition, the induced velocity 
within the body must be zero. Thus, in contrast with the distributions 
previously considered, the Poincare transformation offers the possibility of 
matching the boundary conditions on both the interior and exterior bound- 
aries of BLW. 

Let us suppose that vorticity W = curl v is present in a domain D, 
bounded by a closed surface S, and denote the domain exterior to D by E. 
Here v denotes the velocity vector of the fluid flow. We shall distinguish 
between a fixed point P(x,y,z) at which induced velocities are calculated, 
and variable points of integration, Q(&,n,7t). The position vector from P 


Eq O 2S r and has the magnitude r 


PQ PO) 


The velocity induced by a vorticity distribution can be expressed 


either by means of the Biot-Savart Law, 


a al WT po 
Vp = Dan ; 3 dt (51) 
DEO 
or in terms of the vector potnetial, 
aa ey ese ea (52) 
P 47 P ie 
D ao 


20 


We shall need to distinguish between the vector operators 


= ©) 


da C) ae) 
Ve = Ox 


= == = 
ud ae vk ap ad WS be 


E 7 Ps 5 
The Poincare transformation is as follows: 


ee 1 am 1D) 
+ 4 277 Vp» Pons (53) 


where n is the unit vector in the direction of the outward normal to S. 


First suppose that P is in E. Then we have 


<a 
lac} 
x 
Kl 
td |X 
© |S 
Qu 
q 
I 
J 
ine} 
x 
Tm a 
J 
x 
Ce 
5 
rd | 
ce) 
SS 
I 
p<) 
a 
im) 
hd ie 
OC 
Se 
x 
<i 
(ae 
Qu 
4 


= V, * [ Beas + i ae WV, u |e 
f 20 A PQ PQ 

-v,* [ass [ 5-0, dar (54) 
S Ee D EW 


Zl 


But 


_ [ Seer- |= ds G»5)) 
D EQ S FO 


Substituting (55) and (54) now yields (53) for P in E. 


Next consider the case that P is in D, and let ¥ denote the volume 


and So the surface of a small sphere of radius To about P. Then we have 


wpx [Be ar-v,~ f Seacteyx [ Man (56) 
IP r P r P 16 

0 
where D' =D - Vo: The last integral in (6) is proportional to the 


velocity induced by the vorticity within ¥o> according to (52) and hence 
must vanish as the radius of the sphere approaches zero. Also, by 


Equations (54) and (55), we have 


ae (57) 
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and 


oe ee | ae (58) 
PQ PQ S PQ S PQ 


The convention on the positive sense of n requires that it be directed 


inwards on the sphere S For the term in (5/7), we have 


0° 


v. x we ag = Tiei\. @ey ages | Sew? & 
Pp r Pir D. 0 
ge EM S PQ S z 


0 0 0 Q 
=-[ ak a0 + on, (59) - 
So 
2 5 
where dS = > dQ; and in (58), 
V [ Sba-- [ aa. sa (60) 
P ie 
S 2 S 
0 0 
Also we have 
lim a | a2 dt = V, | CY ar (61) 
£970 eee FO 


since the velocity field of a volume distribution of sources is continuous. 


Collecting the results in Equations (56) to (61) now yields (53) for P in D. 
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Finally, when P is a "smooth" point of S, i.e., a point at which the 
tangent to S is continuous, we introduce a small hemisphere about P, of 


radius Tos and apply the Poincare transformation to the so-—diminished 


volume D', and to the bounding surface, consisting of the hemispherical 


surface So and the remainder S' of S. The proof is similar to that for 


P 


Figure 8 


P in D. In Equations (57) and (58), we need only to replace S by S', and 


in Equation (59) Amv, becomes 2Mv_, or, indeed, av. if P is a corner point 


py 
of S of solid angle a. Instead of (61), we need 


IP 


and 


which are verified by observing that the velocity fields of the surface 


distributions of vorticity, n <X v, and of sources, v ° n, are continuous. 


Equivalent Irrotational Flow From Poincare Transformation 


Let us apply the Poincaré transformation (53) to a case in which the 
flow is entirely due to the vorticity in a domain D, bounded by a closed 
surface S. We shall seek to express the velocity at a point P of E as the 


gradient of a velocity potential. 
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The first term of (53) is seen to give the velocity Vir due to vortic— 


ity in the form of the vector potential (52), as ATV 5 Since the assumed 
flow is solenoidal, the third term of (53) is zero. We then have, from 


G3); 


Amv» = Vp x | — ds - V, i = ds (62) 
su ae 


This expresses the velocity in terms of that induced by a source distribu— 
tion of strength v_/ 4m and by a vortex sheet of strength avs Wel both 


on S. Here we have expressed the velocity vector v = nv, + SV 5» where s 


is the unit vector in the direction of the projection of v on the tangent 
plant at Q. Then, putting n X s = t, we obtain the form given above. 

In order to express the field of the vortex sheet as the gradient of 
a potential, let us define a function > harmonic in D, which on S 


satisfies 


im x vy = < Ws vy = Vo, (63) 
We have then 
3b) 
AEE ate eee S (64) 


and hence, by integration along s, the values of > on S may be presumed to 
be given in terms of the known values of Veo Thus the boundary condition 
(64) sets up a Dirichlet problem for determining 1° 


Since V X v, = 0, application of the Poincare transformation to Vv 


1 P 
taking (63) into account, yields 
- = v°n 
v_ x | a ds = V, i = ds (65) 
Pe ene tO f “EQ 
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Hence (62) becomes 


n° (v_-v) 
Vv -—V, a (66) 
S 


a velocity field due to a source distribution of strength 
GS i © (y/o (67) 


The form of Oo shows that n ° Vy is the normal component of the velocity on 
the interior side of the source distribution on S. Uniqueness of solutions 


of Neumann problems on S then shows that v, is the irrotational velocity 


1 
field in D associated with the source distribution © on 5S. 

An alternative source distribution on S can also be found directly as 
the solution of the exterior Neumann problem for the given values of ven 
on S, as was done previously in considering the displacement effects of the 


boundary layer. Applied to the vorticity field BLW in the flow about a 


body, the present approach requires that a composite bounding surface, 


=S _+T+ 
SS 8,0 Ut A 


be used, where Sp is the surface of the given body, and A is the transverse 
surface of the truncated wake. If the transverse section is taken suffi- 
ciently far downstream, the effects of the source distribution on A may be 
neglected, and the Neumann problem could be formulated as a pair of simul- 
taneous Fredholm integral equations of the second kind. The resulting 
source distribution, however, would not coincide with that given in (67). . 


In the present treatment, the value v * n= 0 nS. is preserved, whereas, 


B 
in deriving (67), the tangential component v xX n = O on Spe in accordance 
with the nonslip condition, was preserved. In the former case, 53 remains 


a stream surface, in the latter, an equipotential. 
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Betz Method for Determining Viscous Drag 

Betz [5] and Landweber and Wu [6] use equivalent irrotational flows to 
derive formulas for the viscous drag of a body in terms of measured values 
of pressure and velocity at a transverse section of the wake. A refine- 
ment of these derivations, in which additional wake characteristics are 
taken into account, will now be presented. 

The body is taken at the center of a circular channel of large radius, 
and is at rest in a uniform stream of velocity U in the positive x- 
direction. The disturbance velocity components in a rectangular (x,y,z) 


coordinate system are (u,v,w), and p denotes the pressure. 


Figure 9 


We select a control surface consisting of the transverse sections AB, 
far ahead of the body, CD or S a moderate distance behind it, and the 
portion of the channel wall lying between these sections. On the section 
AB, designated S 
(U 


0? the pressure is the constant Po> and the velocity is 


0? 0, 0). Application of the momentum theorem to this control surface 


yields the expression for the body drag D, 


D = \ {py-P-P [ Uytu) Uy) ds (68) 
S 


in which 9 is the mass density of the fluid. If the wake is turbulent, 


Reynolds stresses will be present, but these can be taken into account most 
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efficiently by averaging the resulting expression derived for the drag. 


In terms of the total heads, defined by 


e dL 2 < de DoD. ee? 
where g is the acceleration of gravity, (68) becomes 
od a 2 2P 22362. 
D = | og (H)-H) - 7 [ (U,+u) Up -v -w |} ds (70) 


S 


We now consider an equivalent irrotational velocity field (Uy + Uy >Vy> 
W1)> with pressure Py generated by a volume distribution of sources of 
strength HW in BLW, such that (uj 5V, >¥,) = (u,v,w) on T, the outer boundary 
of BLW. We again apply the momentum theorem, to the flow within the same 
control surface generated by this distribution of so-called Betz sources, 


to obtain the expression for the force on the sources within the control 


volume, 


te DPD Ie Lao) 
D =- 7 { [(U,tu,) Up Wig 0h ] ds (71) 


the term corresponding to the difference in heads vanishing since the 
field is irrotational. 
Another expression for the force on the Betz sources is given by the 


Lagally formula 


D. = -— 4Tp { u(Uptu, ) dt (72) 
ye 
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where * denotes the control volume. This is not the usual application of 
the Lagally theorem, which gives the force on a closed body. Here it 
represents the reaction on the Betz sources due to the flux of momentum of 
the source-generated discharge. 

Since the flux through the bounding surface T is the same with Betz 
sources as in the solenoidal (divergenceless) real flow, the difference in 
flux for the flows through the area A of S intersected by the wake is at- 


tributable to the Betz sources. The Gauss flux theorem then yields the 


At | LU dt = f oeo dS (73) 
+ A 


Far downstream, the section S will be denoted by S, and the wake area 


formula 


by A.. In terms of the flux Q across this wake area, we have the well- 


known formula for the drag, 


D = 6 WO q-- f was (74) 
A 


[oo) 


A similar expression for the force on the Betz sources, obtained from the 


asymptotic form of (7/1), is 


But 

J uy dS = | (u,-u) dS + i] u dS 

S08 So So 
and since uy = u except within the wake and, by continuity, i u dS = O, 
we obtain S 


co 
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and hence, from the asymptotic form of (73), 


De = - 410 Up | u dT 
BLW 


Comparison with the asymptotic form of (72) now yields 


4m Wale = © 


This is derived by neglecting uy 


Betz sources are concentrated near the body, so that u 


in comparison with u in (6), since the 


(75) 


(76) 


(77) 


(78) 


(79) 


1 diminishes as the 


inverse square of the distance, while, for a turbulent wake, |u| diminishes 


as the inverse 2/3-power of the distance. 


ca ae . Ege 
Put ¥+ ¥ Vonw and define a mean value uy of uy 
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(80) 


We then obtain, from (77), 


But, by (73) and (78), 


J uae | par- | ear= 2-2 [ cw as 


y' BLW ¥ A 


Then, by (74), 


2 ge | ee av eyes (81) 
a 4t pU5 les 
¥ A 


Hence, by (73) and (81), (72) becomes 


u 
- Sh 
Dy =-0 i (u,-u) (Ujtu,) dS + Uo D (82) 
A 


A formula for the viscous drag D can now be obtained by subtracting 
the expression for dD, in (71) from that for D in (70) and then substituting 


for Ds from (82). Observing that the resulting integrand is nonzero only 


over the wake area A, we obtain the result 


0/2 


- De? 2 2 
[ Pee a0te po oye aye AW V5 8, ] ds (83) 
1-4 
0 


Sill 


In applying (83), total head tubes which measure a total head Hy 
given by 


pg H = p+ 5p [(Uytu) AQ], O<A<1 (84) 


where A is a calibration constant, (see Reference [9]), and velocities 


1/2 
a= | 2 cup) | - U, (85) 


are used for H and u. In terms of HL and ui. with the small difference 


between u,, and u neglected in higher-order terms, (83) becomes 


p = oe { [2g(H)-H )+(u,-u_) (u,tu-2u,) 


1-u,/U, \ 


wv Paw tA (wv?) I dS (86) 


When the wake is turbulent, the mean value of D, obtained by replacing 
(u,v,w) by (utu',vtv',wtw') and averaging, where (u',v',w') denote the com- 


ponents of the turbulent velocity fluctuations, becomes 


- 2 2 
D = —— i [2e(Hy-H,)+(a,-u,.) (atu, -2 )vy ai 


A 


+h Cae) | (Panty | ds (87) 


where Ie and u_, now refer to their mean values. The Reynolds stress terms 


combine into (ijn for isotropic turbulence, and would hence be 
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negligible for 4 = 1/2. Jin Wu's measurements [7] indicate that the tur- 
bulent stresses would contribute about 2 percent to the calculated drag 
with } = 0; but his actual value was i = 0.5, with which the turbulence 
terms became negligible. We shall assume that the turbulence stresses in 


(87) may be neglected. We shall also neglect the terms - v ave W z +X 


il alt 
(ae) since these are small and partly cancel each other. 
Still unknown are the terms uy and a in (87). For estimating Uy> 
we shall assume that uy depends upon z alone, and is given by u, = u, (2), 


the measured value of u at the wake boundary T. According to the defini- 


tion of u, as a mean of u, in ¥', we observe that the mean is weighted by 


i 1 
the value of the source strength U which, together with Uy> diminishes to 
zero aS x > ©, This suggests that the values uy = 0 and wy = uy at A 
can be used to obtain bounds for the drag formula, the "true" value lying 
closer to the bound given by u, = u.,. Because u, occurs both in the 


iL K 1 
integrand and in the denominator of the expression for the drag, it is not 


immediately evident which of the two bounds is the larger. Denoting these 
bounds by D, when i = © ad WD, when vw, = une and applying the afore- 


1 iL 2 il 
mentioned approximations, we obtain from (87), 


_ @ ee ZB 
D, BEG { [2g (Hp Http un ] ds (88) 
A 
and, with up F ws 2u, replaced by u,7 u, in the last term, 
D5 = = fA [2¢(H,-H_)-(u_-u 7] ds (89) 
om fU O m E m 
is) A 


The latter form was given by Tzou and Landweber [8]. 

That dD, € Dy is indicated by the following argument. Since the source 
strength ) represents the displacement effect of the boundary layer and 
wake, it is a positive quantity of total strength given by (78). Conse- 


quently, according to (77), u, cannot be of one sign throughout BLW. Over 


il 
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most of the boundary-layer edge of T, up > 0. In the region of the rear of 
the body and the near wake, there would be some pressure recovery (especial- 
ly if separation has not occurred) and u, on T would there become negative, 
and gradually approach zero with increasing downstream distance. This 


suggests that u, in (89) and uy in (86) are negative. Since the difference 


E 
in heads, Ho - Ho contributes about 90 percent of the magnitude of the 


integrand in (88) and (89), and a, << U (89) may be written as 


0? 


= 2 5 i a) paste 
dD, 35 (+ 7 ) I 28 (H HD ds 5) { (up u) ds 
0 A A 


which yields 


oe ree ® _ df 
D, - DL SiG { g (Hy HD ds fo) [ u, (up uD ds 
v A A 


u 
= 9 r || [g(H)-H_,)-Up(u,-u)] ds (90) 
A 


Applying the expression for Hn in (85) and 
cS i 2 
Pg Hy = Py t > PCUptu,) 


in (90) and neglecting the term u aa oo we obtain 


ica] 


(Pp-P_,) dS > 0 (91) 
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since ar <a ee Amand accondinom ton ammWin smd atzais[slOy| = Pe is a maximum at 


the center of the wake, so that Diag: 2 0. The data given in Reference 


p 

E 
[8] indicate, however, that, at a section at 0.6 of the length behind the 
stern of a ship model, Up is zero within the accuracy of the measurements. 
This indicates that the simpler expression (88) is suitable for computing 


the viscous drag from wake-survey data. 
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